c-axis transport in highly anisotropic metals: role of small polarons 
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We show in a simple model of interlayer hopping of single electrons, that transport along the 
weakly coupled c-axis of quasi-two-dimensional metals does not always probe only the in-plane 
electron properties. In our model where there is a strong coupling between electrons and a bosonic 
mode that propagates in the c direction only, we find a broad maximum in the c-axis resistivity at 
a temperature near the characteristic energy of the bosonic mode, while no corresponding feature 
appears in the ab plane transport. At temperatures far from this bosonic energy scale, the c-axis 
resistivity does track the in-plane electron scattering rate. We demonstrate a reasonable fit of our 
theory to the apparent metallic to non-metallic crossover in the c-axis resistivity of the layered 
ruthenate Sr2Ru04. 
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Many metallic systems of current interest have highly 
anisotropic electronic properties with one or more direc- 
tions much more weakly conducting than the others. Ex- 
amples include the cuprate metals*, some ruthenates*, 
and many organics 3 -, where a measure of this anisotropy 
is the resistivity ratio, p c /p a b, which can be of order 10 3 
to 10 5 . (Henceforth, c will denote the weakly coupled di- 
rection with a and b the more conducting.) Many of these 
materials have unusual metallic or superconducting prop- 
erties which have been attributed to strong correlation 
effects in the ab plane. Although the theoretical under- 
standing of strongly correlated electrons is incomplete, 
the electron anisotropy can be used advantageously to 
give insights into the in-plane physics. By assuming that 
the coupling up the c-axis occurs only by weak single- 
electron tunneling, many workers^fi have shown that 
c-axis resistivity depends only on the in-plane electron 
Green's function. 

This insight makes all the more puzzling the obser- 
vation that in many of these materials the resistivity- 
anisotropy ratio is strongly temperature dependent with 
a cross-over from metallic behavior to insulating be- 
havior in the out-of-plane direction only. For example, 
Si'2Ru04 2 shows T 2 resistivity both in the plane and out- 
of-plane for T < 20K but p c shows a broad maximum at 
around 130K and above this temperature the resistivity 
starts to decrease with increasing temperature (see Fig- E] 
below). In marked contrast, the in-plane resistivity shows 
dp/dT > for all temperatures. Evidently the low tem- 
perature metal is a 3D Fermi liquid, so this begs the 
question how the out-of-plane transport is being blocked 
while the in-plane transport remains metallic. In this pa- 
per we consider how anisotropic coupling of the in-plane 
electrons to bosonic degrees of freedom can give rise to 
this effect. We show that if the electron-boson coupling 
is stronger in the inter-plane direction then, at the char- 
acteristic energy scale of the boson mode, there can be a 
maximum in the c-axis resistivity with no corresponding 
transport feature for currents in the plane. 

Other approaches to this puzzle have been offered pre- 
viously. One possibility is that the in-plane metal is a 



non- Fermi liquid^ or otherwise reflects unusual in-plane 
physics^ such as superconducting fluctuations^. These 
approaches have generally been argued to apply to the 
cuprate metals, but our motivation comes from the lay- 
ered strontium ruthenate compounds where the appeal 
to unusual in-plane physics is less well grounded. Alter- 
natively the c-axis coupling may consist of more than sin- 
gle electron tunneling: with disorder and boson-assisted 
hopping 11 or perhaps inter-plane charge fluctuations^. 

In this paper we also go beyond single particle hopping 
as the only coupling between the planes and consider elec- 
trons coupled to a bosonic mode. However, in contrast 
to the above work, we use a canonical transformation 13 
to treat the electron-boson coupling, and thus can deal 
with the strong coupling regime. Our assumptions are, 
first, that there is a well defined separation of energy 
scales: t a b » uo > t c , where t a h is the in-plane ki- 
netic energy scale, wq is the characteristic energy of the 
bosonic mode {e.g. the Debye temperature) and t c is the 
single particle tunneling between layers. This assump- 
tion is certainly well founded for Si'2Ru04, as shown, for 
example, in Ref. [l4|. With this assumption it is suf- 
ficient to compute conductivity to leading order in t c . 
Our second assumption is that electron-electron correla- 
tions dominate over the in-plane coupling to the bosonic 
mode such that the formation of a small polaron in the 
planes is suppressed by the electronic interactions. We 
will therefore be considering a limit where the quasipar- 
ticle excitations of the system are insensitive to bosonic 
modes whose wavelength in the plane is short compared 
to the size of a quasiparticle (assumed to be large) . With 
these assumptions we show that anisotropic coupling can 
lead to differences between the in-plane and out-of-plane 
conductivities that are consistent with the experiments. 
A broad maximum can appear in the c-axis resistivity p c 
as a function of temperature, while no such feature exists 
in the in-plane resistivity p a b- Only away from the broad 
maximum, does p c track the intrinsic in-plane scattering 
rate. We shall also show that our model calculations can 
fit reasonably the Sr 2 Ru04 data of Tyler et al£, using 
physically plausible parameters. 
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Strictly our assumption that the electrons couple 
strongly only to out-of-plane bosonic modes should be 
regarded as phenomenological, and in this paper we com- 
pute the consequences of such an approach and compare 
with experiment. Nevertheless we can offer some justifi- 
cation for it. In a quasi-2D metal, the anisotropy of the 
crystal structure has important consequences for both the 
electronic structure and the bosonic modes the electrons 
can couple to. Usually the anisotropy in the characteris- 
tic energies of phonon modes are much smaller than those 
than the electronic structure: the ruthenate that we will 
be concerned with is thought to be an ionic (not covalent- 
bonded) crystal. Where significant anisotropy can arise 
is in the strength of the electron-phonon couplings for 
these different phonons. Krakaueref. al*^ have found 
in an LDA calculation that the electron phonon cou- 
pling is much stronger when atomic displacements are 
perpendicular to the Cu-0 plane in LSCO, while Kim 
et. aZ.— and Grilli and Castellaniii found that strong 
in-plane electronic correlation suppresses in-plane elec- 
tron phonon coupling. Usually strong electron-phonon 
coupling leads to the well-studied small polaron regime, 
where electronic motion becomes affected by the accom- 
panying large cloud of phononsi2*i£. (For a review, see 
flflj.) In this paper, our assumptions will lead to a small 
polaron forming only between the layers. 

Our model Hamiltonian is then: 

H = J2 Hin) + H c + H c - b + H h , 

n 

H c = f d 2 xY,tcci +1 (x)c n (x)+U.c. . (1) 

J n 

The model system consists of a stack of 2D planes de- 
scribed by H^ 1 ' containing arbitrary intra-planar inter- 
actions. In this work it will be sufficient to characterize 
these planes by a phenomenological 2D Green's func- 
tion. In H c , electrons are assumed to hop from one 
plane directly to a neighboring plane only, at the same 
in-plane coordinate x. For simplicity, we have omitted 
the spin index on the electron operators c n {x). Many 
approaches*^ assume only this coupling between the 
planes in a quasi-2D system. We go beyond that by in- 
cluding a bosonic mode that also couples the layers 

H e -b = ^ ^ M q exp(iqnc)p n (a q + oO . (2) 

n q 

Our analysis can be applied to any gapped neutral 
bosonic mode with long wavelength in the plane, such 
as a magnetic collective mode, but we have in mind 
the (strong) coupling to a phonon mode with displace- 
ment in the c direction and long (infinite) wavelength 
in the plane. Here q is along the c-axis, with c the 
inter-layer distance. In this limit, only the total charge 
p n = J d 2 x(y n {x)c n (x) of the n'th plane couples to the bo- 
son creation and annihilation operators a q , a q so, by con- 
struction, these bosons do not dominate charge transport 
within the plane. 



Even with this limiting electron-phonon coupling, 
interlayer hopping can affect in-plane correlations by 
spreading the in-plane correlation from one plane to an- 
other, especially since the in-plane correlations are sup- 
posed to be strong for a pure single layer in the cuprates 
and also the ruthenates. However, this requires the inter- 
layer hopping to be coherent, which means the temper- 
ature must be low compared to t c : otherwise, the c-axis 
transport is diffusive and cannot propagate correlations 
between the layers. In this paper, we have assumed a 
phenomenological form for the in-plane spectral weight, 
and so we cannot address the issue of the strongly corre- 
lated three dimensional state at T -C t c . 

The Hamiltonian for the c-axis bosons, Hh = 
*52 q u q a q a q) defines a dispersion u q . In this paper we 
shall look closely at the Einstein phonon uj q = loq. This 
is both because optical phonons tend to have little dis- 
persion relative to acoustic ones, and also we are able to 
obtain analytic results for this case. We shall also dis- 
cuss, somewhat more qualitatively, the opposite case of 
a generic dispersion where the phonon density of states 
does not have any sharp features. For the data fitting to 
Sr2Ru04, we will use a simple form that allows us to test 
if the putative mode does indeed disperse. 

Of course electron-phonon systems have been well 
studied, but our findings will be shown to be clearly dif- 
ferent to ones found in the small polaron models stud- 
ied in the classic papersiiiaiii: small anisotropy in 
electron-boson coupling cannot lead to the properties we 
shall demonstrate below. In particular, at temperatures 
greater than the temperature of the maximum in p c (T), 
our model has p c (T) tracking the intrinsic in-plane scat- 
tering rate while in standard small polaron theory (even 
if generalized to some degree of anisotropic coupling), 
p c (T) keeps on decreasing (exponentially) with increasing 
T . In short, in our model, the small polaron broad max- 
imum in p c (T) is "grafted" on top of the background in- 
plane scattering. We note that KornilovitcbSS has stud- 
ied a similar anisotropic model: he did not calculate the 
d.c. conductivity, and, with the effective mass approach 
employed, he cannot access the non-metallic regime [be- 
yond the broad maximum of p c {T)\. After completion 
of our theory, we became aware of the related work of 
Lundin and McKenzia^i who also studied a small polaron 
model for anisotropic metals. Their model is different to 
ours in an essential way: they take the boson mode to be 
uncorrelated from one layer to another, while we have a 
boson mode that propagate coherently along the c-axis. 

The key physics we are considering is the effect of 
strong electron-phonon coupling on the charge transport 
in the weakest direction, the c-axis. Physically, the mo- 
tion of the electron is accompanied by the emission and 
absorption of a large number of phonons due to the 
strong coupling, forming the so-called small polaroid. 
We treat -ff e -b exactly by the canonical transformation, 
H = exp(—S)Hexp(S), a straightforward generalization 
of the transformation well known in the small polaron 
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,13,19 



S = - f d 2 xJ2 ^e^ nc ct (x)c n (x) (a q - al q ) . (3) 

n,q 



Then, H = Y, n H (n) +H c + H h , and 



s - = -E 



M, 



,f |2 



-Pn 



53 J d 2 x tJ n+1 (x)c n {x)xl +1 X n +YL.c. , (4) 



where 



(5) 



Thus strong electron-phonon coupling leads to (1) a 
renormalization of the in-plane chemical potential (first 
term of S c ), which we shall henceforth ignore, and (2) an 
effective vertex correction for the c-axis hopping t c (sec- 
ond term of H c ). At temperatures lower than the charac- 
teristic energy of the boson, electrons can hop coherently 
from one plane to another, with some suppression due to 
the "dragging" of the boson cloud. With increasing T, 
the hopping electron gets inelastically scattered by more 
and more bosons: inter-layer hopping is now more diffu- 
sive. Hence a crossover to non-metallic behavior occurs. 

We consider the d.c. conductivity and the orbital 
magnetoconductivity for fields in the ab plane. For 
the conductivity, since the charge in the n'th plane is 
Q c n = e J d 2 xcl l (x)c n (x), the current in the c-direction is 
just jn — ®tQn = — i[Qn'H]- After the canonical trans- 
formation, jn — ~ i[QniHc], and so: 

tin) = *e*c J d 2 x ( [ci +1 (x)c n (x)xl +1 X n - h.c] ) . (6) 

The orbital effect of a magnetic field can be included 
using a Peierls substitution 2 ^. The mathematical details 
of our derivation follow our analysis reported in Ref. |2j 
Using the Kubo formula and expanding to 0(t 2 ) gives 
the conductivity, i.e., the linear response to an applied 
electric field in the c-direction a c = j c /E c : 



= lim£ 7T 



E 



d kn / dr 



xe^U(r)($ D )$ h T)G%3fa +&,r) (7) 



where C/(T) = (4 +1 (r)X„(r)Xt(0)X n+1 (0)) , (8) 

Gn D \k\\) is the in-plane (dressed) Green's function for 
plane n at in-plane momentum fcii , and qs = ec x B 
with c being the inter-plane lattice vector. The temper- 
ature dependent factor U (r) in Eq. [7| leads to the c-axis 



conductivity being no longer simply a convolution of the 
in-plane electron Green's functions and we explore the 
consequences in the following. 

To make further progress we use a phenomenologi- 
cal in-plane electron Green's function in the spectral 



and 



representation: G n '(k\\,u) m ) = f dz iulm _ z , 

U(y n ) = J dz M^zr , then the Matsubara sums can be 
done. We shall take the in-plane electron spectral func- 



tion to have the form: A( 2D '(fc, u) — — ^_ e ^ 2+r ^ T ^y i 

where we assume that the scattering rate T(T, to) has no 
in-plane momentum dependence. Then, integrating over 
fen , the zero-frequency, zero-momentum c-axis conductiv- 
ity becomes: 

a c (T, B) = ^Re f dudv^ ~ n *< W + v) 



ITcT 



D(u,T) 



[(v F q B ) 2 -{v + ir(T, w) + iT(T, uj + v)f] 1/2 



(9) 



where D(u>,T) — uib(ui)nB{uj) [1 + ns(w)], Uf,b are the 
Fermi and Bose distributions and N is the electron den- 
sity. D(ui,T) can be calculated 2 ^ following Reffl9l. We 
now analyze this equation in detail for Einstein modes 
and, more qualitatively, for a general dispersing mode. 

For Einstein modes, u> q = uq, the function D(u>,T) 
is made up of delta functions at the harmonics u> = 
nujQ, n = 0, ±1, ±2, . . It can be shown that for exper- 
imental temperatures, one needs only the first few har- 
monics. We now discuss the three temperature regimes 
(T <C ujq, T < u>q, T ujq) separately, and show that the 
low and high temperature regimes basically tracks the 
in-plane scattering rate, just like in previous studies of 
interlayer transport 4,5,6 . Only in the regime where the 
temperature is near the boson scale uiq can the electron- 
boson scattering dominate over the in-plane scattering 
contribution, and instead, leads to a small polaron like 
broad maximum in resistivity. 

At low temperatures (T <C wo), the asymptotic form 
for the c-axis resistivity is the same as for band elec- 
trons, but with an effective hopping parameter t^. s = 
t c e~ A2/2uJ o, where A 2 = 2 ? 2|M g | 2 (l - cosg) character- 
izes the strength of the electron-phonon interaction: 



e 2 Nt 2 c 
2c7rr(T,0) 



e-^M 2 (T«w ). 



(10) 



For T ^> uj , then a c ~ T _7)_ 2, using the scattering 
rate T(T) oc T v . (rj = 2 for a Fermi liquid, while for the 
marginal Fermi liquid, r) = 1.) Thus the low temperature 
-6 = conductivity directly probes the in-plane scatter- 
ing, while the high temperature conductivity still reflects 
the in-plane scattering with an extra T -1 / 2 factoid. 

However, beyond a critical phonon coupling (A c ) there 
is an intermediate temperature region where a broad 
maximum appears in the c-axis resistivity at tempera- 
ture T max (Fig. [TJ. Beyond T max , the resistivity dips to 
a (broad) minimum at T m j n before finally joining onto 
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the asymptotic T ^> uj regime mentioned above . We 
estimate, by using only the zeroth harmonics n = 0, that 
if the in-plane scattering rate has the form T(T) = aT n , 
then the maximum appears when 



A>A c «^ 0V V0-4. 



(11) 



The position of the resistivity peak for A > A c is found 
to be T cr it ~ luq/3.0 . For A only slightly larger than A c , 
is given by: 



u 



2Ti 



max 2X' cr it 
min 



± 



\ 



JL 

0.2 



up 
A,. 



(12) 



Note that both T max and A c depend mainly on phonon 
parameters; one can show that the in-plane scattering 
rate enters only in the form of the exponent rj. Hence the 
width of the broad maximum is governed by the scale ojq. 




T/co, 



FIG. 1: The zero-field c-axis resistivity with in-plane marginal 
Fermi liquid scattering F = u> coth (w/T) for electrons coupled 
to c-axis Einstein phonons for different values of electron- 
phonon coupling (A/loq). Inset: as above but for a dispersing 
mode Ug = u>o[l — 0.4cos(g)]. 

For a more general dispersing bosonic mode, D(co,T) 
can be calculated approximately^ assuming strong 
electron-phonon interactions. Using the usual approxi- 
mation [np{v) — np{y + cj)]/cj ~ —dnp/dv ~ S(v) valid 
at low T when T(T, u>) varies slowly, we find 



Cc(T,B) = 



dv 



xRc 



<<<'T J y/irffT) sinh fa 

CX P ( — 47= (T) ) 



[{v F q B ) 2 - (v + 2iT{T, v)) 



21 1/2 



(13) 



where 



, s v-^ 2 sinh 2 {%) „ 



sinh 



2T 



and F q = 2(1 - cosg). 

For illustration, we consider uo q — oj [1 — 0.4cos(g)], 
and set the electron-phonon coupling M q to be q- 
independent. Eg. 1131 is then evaluated numerically and 
plotted in the inset to Fig. ^ A = ~Yli q F q u! q /2 charac- 
terizes the strength of the electron-phonon interaction in 
this case. The result is qualitatively similar to that of a 
non-dispersing mode. 

Now we consider the magnetic field dependence of the 
conductivity for the Einstein mode. At low tempera- 
tures T < wo we find that the magnetoresistance reflects 
the usual cross-over from quadratic to linear field de- 
pendence^ at a scale determined only by the in-plane 
electron scattering rate 



Ap c = [ Pc {B) - p c (0)]/Pc(0) = ^l + [a; c /2r(r,0)] 2 , 

(14) 

where lj c = vpecB is the cyclotron frequency. At higher 
temperatures (T > loq), the weak-field magnetoresistance 
is quadratic in field Ap c oc (w c /r c ff ) 2 , and defines a new 
scattering rate, r e ff, depending on both the electronic 
scattering and the phonon frequency. Taking into ac- 
count only up to n = 1 harmonics, with To = T(T, 0), 
and Ti = T(T, ujq), we find 



1 off 



1 

2r 



4r 



4Ti 



1 2r 1 (4r;-3o;g) 



(15) 



2sinh(w g /2T) ' 



As for the zero field resistivity, the magnetoresistance for 
the dispersing mode case is qualitatively similar to that 
for the Einstein mode and is positive for all temperatures. 

We now consider in more detail how this model might 
apply to real materials and, in particular, the ruthen- 
ate systems. The 2D electronic nature of these materials 
reflects the crystal structure and so it is natural that 
the electron-phonon interaction should be anisotropic. 
While there are to date no calculation for the electron- 
phonon interaction parameters for the ruthenates, in the 
iso-structural LSCO cuprate family, Krakauer et. alii 
have calculated that there is a strong electron-phonon 
coupling only for modes corresponding to atomic dis- 
placements perpendicular to the Cu-0 plane, (partly) 
because of weak screening of the resulting electric fields 
in this direction. Moreover it has been argued that in the 
perovskite structure the coupling to c-axis vibrations is 
enhanced^. There exist optical phonons with the appro- 
priate symmetry for c-axis transport 28 in Sr2Ru04, and 
experimentally the broad maximum in the c-axis resistiv- 
ity has been linked to a structural phase transition 2 ^ in 
Cai.7Sro.3Ru04 at around the broad maximum tempera- 
ture. Thus one should consider the possibility of electron- 
phonon interaction affecting the c-axis transport. In par- 
ticular, both the Sr2RuC>4 and Cai.7Sro.3Ru04 systems 
cxhibitSiSL qualitatively this broad maximum structure 
found in our simple model, in the c-axis resistivity near 
to their characteristic (c-axis) phonon energy. 

A full quantitative modeling of Sr2Ru04 is beyond the 
scope of the simplified model presented here. Despite all 
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FIG. 2: Fitting the zero-field c-axis resistivity to the data of 
Tyler et.al.—. See the text for the explanation of the theoret- 
ical curve (thin line) and the parameters used for the fitting. 
The data are plotted as points: for clarity, we have plotted 
only 3% of the total number of available data points. 



these simplifying assumptions, our model does a reason- 
able fit to the resistivity data of Sr 2 Ru04 as we shall 
show, perhaps indicating some degree of universality of 
the mechanism for c-axis transport studied in this paper. 
Here for reference, we list the simplifications we have as- 
sumed in our modeling: 1) the in-plane scattering rate 
can be deduced directly from p a },, 2) the main interlayer 
coupling is single particle hopping from one layer to an- 
other and 3) a simple phonon dispersion is employed (see 
later). 

For 1), this amounts to ignoring vertex corrections to 
in-plane transport. In Fermi liquids, the vertex cor- 
rection leads to an extra cosine of the angle between 
in- and out-going electrons thereby correctly penalizing 
back-scattering in conductivity. However, the qualita- 
tive trend is still correct. As mentioned already, we have 
taken an in-plane spectral function that has no in-plane 
momentum dependence. This ignores the complicated 
multiple Fermi surfaces in Sr2Ru04. For 2), this means 
we ignore interlayer coulomb interaction. Also, we ig- 
nored the multiband nature of Sr2Ru04 and the depen- 
dence of t± on in-plane momenta. For 3), we have taken 
the limiting case of an optical phonon mode where the 
atomic displacements are perpendicular to the plane, and 
all the atoms in the plane move together. We envis- 
age that just as in LSCO, there will be strong electron- 
phonon coupling only for phonons propagating mainly 
in the c direction- 5 ". Now in reality, because of the non- 
trivial perovskite structure of Sr2Ru04, these modes that 
can affect c-axis transport will have relative atomic dis- 
placements both in-plane and out of plane. But this 
makes electron hopping from one to another plane even 
more difficult, as the polaron has to create disturbances 
both in-plane and out-of-plane. 

To some extent, these simplifying assumptions may 



only lead to some quantitative changes in the fit param- 
eters, because the calculation of p c (Eq. 0) involves an 
integral over in-plane momenta and thus averages out 
such dependences. 

We now show in Fig. a fit to the Sr2Ru04 data of 
Tyler et al£. The theoretical curve is the thin continu- 
ous line, the data are the points. The theoretical curve is 
generated as follows: the in-plane scattering rate is ap- 
proximated as being proportional to the in-plane p a \, (also 
taken from the data of [Q): T(T,lj) m A(p ab (T) + p ), 
where A and po are fitting parameters. We take a sim- 
ple optical boson dispersion uj(q) = u>o (1 — Acosg) with 
the boson bandwidth A and the characteristic energy 
ujq as fitting parameters. T(T,u>) and w(q) are then 
fed into Eq. ^3 where the overall scale of p c is found 
by fitting p™ ax ~ 33mi7 cm to the peak of the theory 
curve. Despite the simplicity of the model and the lack 
of knowledge of the exact boson dispersion form, the 
theory curve is almost indistinguishable from the data 
points, except at high Tml The fitting is relatively insen- 
sitive to the precise values of the parameters: some fine 
tuning (mainly, A and A) is needed to get the positions 
of the maximum, minimum and final upturn correctly. 
Furthermore, the parameters used are physically plausi- 
ble: A w 0.67, w « 515K, A/oj « 1.4,A/w ~ 0.1(mfi 
cm) ,po ~ 0.17 mQ, cm. The value of loq is consistent 
with the available phonon data of Braden et. alm&, while 
the size of the boson bandwidth suggests that the bosons 
do disperse in the interplane direction. Note that the 
dispersion u>(q) always enter the conductivity expression 
under a g-integral (see Eci ll3|) . and is the reason why 
our simplified boson dispersion can still model the broad 
maximum in the c-axis resistivity successfully. The size 
of the electron-boson coupling parameter A indicate that 
Sr2Ru04 is just slightly above threshold for obtaining the 
maximum in p c . 

The c-axis magnetoresistance seen in Sr2Ru04 is un- 
usual, becoming negative for both transverse (B||a6) 
and longitudinal (B||c) fields above 80K and maximally 
so at around 120K — coinciding with the peak position. 
In this paper we have only considered orbital magne- 
toresistance and, as might be expected, found it to be 
positive 31 . This is not inconsistent with the data since 
the transverse magnetoresistance is always less negative 
than the longitudinal one. If the change of sign is linked 
to the origin of the resistivity maximum then, within the 
framework presented here, it indicates that the frequency 
of, or coupling to, the bosonic mode is field dependent. 
At present, we have no microscopic picture of how this 
might occur. 

In conclusion, we have shown that c-axis transport in a 
quasi-2D metal does not always probe only the in-plane 
electron properties. Strong coupling between electrons 
and a bosonic mode polarized in the c-direction in a 
highly anisotropic metal can lead to a broad maximum — 
an apparent metallic to non-metallic crossover — in the c- 
axis resistivity with no corresponding feature in the ab 
plane properties. The position in the temperature axis, 
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and the shape of this broad maximum are determined 
mainly by boson parameters and not on the magnitude 
of the in-plane scattering rate. We have discussed the 
potential application of the model to Sr2Ru04 and its 
relative Cai.7Sro.3Ru04. Despite certain simplifying fea- 
tures of the model, the qualitative (and even quantita- 
tive) properties of this crossover in these layered ruthcn- 
ates are captured succinctly. 
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to take this into account within the present framework. 
Other processes may intervene at high temperature, but 
one source of discrepancy at high T is that the experi- 
mental data is the resistivity at constant pressure while 
our calculation is for constant volume. Non-trivial thermal 
expansion at higher T will contaminate the data. 



